Elastic vs. inelastic coherent backscattering of laser light by cold atoms: 

a master equation treatment 
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We give a detailed derivation of the master equation description of the coherent backscattering of 
laser light by cold atoms. In particular, our formalism accounts for the nonperturbative nonlinear 
response of the atoms when the injected intensity saturates the atomic transition. Explicit expres- 
sions are given for total and elastic backscattering intensities in the different polarization channels, 
for the simplest nontrivial multiple scattering scenario of intense laser light multiply scattering from 
two randomly placed atoms. 
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I. INTRODUCTION 



Localization phenomena in disordered systems have 
become a subject of intense research 0, y|, since they 
highlight the fundamental role of interference effects for 
wave propagation. A prominent example is the coher- 
ent backscattering (CBS) of light in dilute, disordered 
media, a hallmark of weak localization. CBS manifests 
itself in the enhanced backscattering of the injected radi- 
ation, due to the constructive interference between time- 
reversed pairs of multiple-scattering amplitudes along a 
given sequence of scatterers, which prevails even after 
disorder averaging. This remarkable effect was observed 
for the first time with light scattering from suspensions 
of polysterene particles |2j, but recently has also been 
reported for laser light scattering from clouds of cold 
atoms 1^ H, Hi ' With such quantum scatterers - possess- 
ing an internal electronic structure which can be probed 
by the scattering light field in a controlled way (e.g., by 
the appropriate choice of the laser frequency and/or of 
the atomic species) - additional decoherence processes 
are brought into play, which fundamentally affect the ra- 
diation transport across the scattering medium. Elastic 
Raman processes on degenerate atomic transitions driven 
by the injected field change its polarization, alike spin- 
flips of electrons scattering from lattice impurities 0- 
Inelastic processes are induced by intense driving of the 
atomic transition, leading to its saturation and a nonlin- 
ear response of the atom ^ , manifest in the emission of 
photons with frequencies different from the one injected. 

Both these decoherence mechanisms reduce the CBS 
intensity, occur in general simultaneously, and are exper- 
imentally very well controlled. Hence, experiments on 
the multiple scattering of coherent radiation on atomic 
scatterers provide an ideal testing ground for the detailed 
analysis of coherent quantum transport in disordered me- 
dia, and of its sensitivity towards various sources of de- 
coherence. Furthermore, if we consider CBS and weak 
localization as a precurser of strong (i.e., Anderson) local- 
ization, decoherence phenomena affecting CBS are likely 



to become detrimental for the latter. Anderson local- 
ization of light, however, is an important experimental 
target, for fundamental as well as for technological rea- 
sons. Consequently, beyond its fundamental interest, a 
detailed theoretical and experimental understanding of 
disorder- and/or decoherence-induced transport phenom- 
ena is highly desirable for possible applications, which 
currently emerge, e.g., in the area of random lasers 0- 

While the impact of elastic spin flip processes on the 
CBS signal is nowadays well-understood, wit h q uanti- 
tative accord between experiment and theory [lOj, non- 
linear processes due to the saturation of atomic transi- 
tions still challenge our theoretical understanding. On 
the one hand, perturbative approaches are - by defini- 
tion - badly suited for the regime of strongly driving 
intensities. On the other hand, exact solutions which 
take into account arbitrarily high multiple-scattering or- 
ders are prohibitive, due to the exponentially increasing 
number of the contributing scattering paths and of the 
coupled internal states of the atomic scatterers. Differ- 
ent approaches are presently persued in the attempt to 
achieve a better understanding of CBS in this parameter 
regime. These range from diagrammatic techniques 
over Langevin equations jl2j | , to a master equation treat- 
ment ^i^, for a small number of atoms. In the present 
paper, we give detailed account of the latter approach. 

We will focus on the scenario set by the first exper- 
imental study of saturation-induced effects on the CBS 
signal, performed with cold Sr atoms 6]. In these ex- 
periments, the injected laser was near-resonant with the 
^5*0 — transition, which has a nondegenerate ground 
state and thus leaves no room for spin-flip processes. 
Consequently, only inelastic scattering could cause deco- 
herence and thus reduce the CBS signal. This was indeed 
experimentally observed already for moderate values of 
the atomic saturation parameter 



(1) 



il is the driving-induced Rabi frequency, 7 half the 
spontaneous decay rate of the excited atomic level, and 
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6 = lul — i^Q the detuning of the injected laser frequency 
LU]^ from the exact atomic transition frequency wq, see 
Fig. [2 

While our formalism to be unfolded hereafter is not 
restricted to the treatment of atomic transitions with 
nondegenerate ground states, this specialization allows 
for a more transparent presentation, and, in particular, 
for a clear identification of the various inelastic processes 
which intervene. 

The paper is organized as follows: The next section 
starts out with a general Hamiltonian formulation of the 
dynamics of N atoms under coherent external driving, 
and coupled to the electromagnetic vacuum. A master 
equation for the time evolution of the atomic degrees of 
freedom constitutes the central building block of the the- 
ory. Explicit expressions for the (back-)scattering inten- 
sities in arbitrary polarization channels are derived, in 
terms of the steady state quantum mechanical expecta- 
tion values of atomic dipoles and of dipole-dipole correla- 
tion functions. Expansion of these to second order in the 
dipole-dipole interaction constant between pairs of atoms 
finally allows us to present analytic expressions for the 
polarization-filtered backscattering signal, assuming that 
double scattering processes provide the dominant contri- 
bution. Accordingly, we restrict our final evaluation to 
the case of light scattering from two, randomly placed 
atoms. Section lllll provides a recipe of how to perform 
the disorder average, before Sec. lIVI presents quantitative 
results for the different polarization channels. Section Ivl 
concludes the paper. 



FIG. 1: Level scheme of a = ^ Je = 1 dipole transition, 
with atomic transition frequency luq, and natural linewidth 
27. The sublevels |1) and |3) have magnetic quantum number 
m — 0. Sublevels |2) and |4) correspond to m = —1 and 
m — 1, respectively. 
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Haf = ft^^ [Kk(i"a)ak,,(D„ -ek.s) 



a = l k,; 



Hal = - 
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II. MASTER EQUATION APPROACH TO 
COHERENT BACKSCATTERING 



Here, Dq-* is the lowering (raising) operator for the 
isotropic dipole transition (see Fig. ^ at the resonance 
frequency ujq of atom a, defined by 



Dq = — e_i(7"2 



- eo(7,3 — e+icr. 



(7) 



A. Full A'^-atom master equation 



The cr^; = \k)^{l\^ mediate transitions between the elec- 
tronic states of atom a, and 



We start with a general formulation of the Hamilto- 
nian describing N identical, motionless atoms with an 
isotropic dipole transition coupled to the quantized pho- 
ton reservoir and driven by a quasiresonant (classical) 
laser field. The total Hamiltonian of the system. 



H = Ha + Hf + Haf + Hau 



(2) 



1 , 

e±i = T-^[ex ± ley), 



eo 



(8) 



are the unit vectors of the spherical basis. In the free 

field Hamiltonian, a'^^ annihilates (creates) a photon in 
the reservoir mode with wavevector k and transverse po- 
larization £k,s, where s is the polarization index. 

The interaction Hamiltonians Haf and Hal are writ- 
ten in rotating wave and dipole approximation. The cou- 
pling constant between atom a, located at point Tq, and 
the vacuum mode (k, s) reads 



contains the free atomic Hamiltonian Ha, the free field 
Hamiltonian Hf , the atom- field coupling Haf , as well as 



2h€aV 



1/2 



(9) 
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where d is a reduced matrix element, eo is the permit- 
tivity of the vacuum, and V is the quantization volume. 
The coupling of atom a to the laser field 

Ei(r) = ELSe'^^^-"-'^^'^ + c.c. (10) 

is characterized by a position-dependent Rabi frequency 

n^^t^e^'^L -. =^le'^^ ''- . (11) 

h 

The figure of merit in our present study is the average 
value of the stationary intensity /(r) with polarization £, 
scattered in a direction close to backscattering — k^: 

/(r)^ lim([£-E(-)(r,t)][£*-E(+)(r,i)]>. (12) 

t — >oc 

E(~/+)(r,t) is the negative/positive frequency compo- 
nent of the source field operators, given by the superposi- 
tion of the retarded fields radiated by all atomic dipoles, 
that is projected onto the polarization vector e, upon 
detection, 

2 ^ 

e* ■ e(+) (r, t) = V e* ■ D„(t„)e-'''--= , (13) 

with ta = t — \v — fal/c, and k the wave vector with 
the wave length of the injected laser radiation, pointing 
in the observation direction (note that all the nontrivial 
spectral information is contained in the time dependence 
of the atomic dipole correlation function) . This expres- 
sion follows immediately from generalizing familiar ex- 
pressions for the far field radiated by a single atom 
to the present case of an atomic cloud, with the cloud's 
diameter much smaller than the distance to the detector. 

The total scattered intensity is then obtained by in- 
serting H13(l and its conjugate into p2(l . and reads, up to 
a prefactor, 

N 

1= J2 ([£•Di][s*•D^])^^e'''•-^ (14) 

a,/3=l 

where 'ss' stands for steady state, and r^f} = — r^. 

B. Coherent backscattering and polarization 
channels 

The structure of Eq. (|14|l . together with (T)), shows that 
all three transitions between each atom's electronic levels 
(see Fig. will in general contribute to the scattered 
light intensity - the relative weights of their contributions 
depend on the observation direction and on the detected 
polarization channel. 

The unit polarization vector e may be chosen either in 
a circular or in a linear basis. For circular polarization, it 
is convenient to use the so-called helicity basis, in which 
the quantization axis is directed along the probe direction 



(a) h\\h (b) lin ± tin (c) h ± h (d) lin \\ lin 

T -6 *f -13 *f 

2 

FIG. 2: (Color online) Elementary configuration for the 
polarization-selective detection of the coherent backscattering 
signal (CBS) (dashed arrows) from laser light (thick arrows) 
scattering on two atoms 1 and 2 (black circles). The dipole- 
dipole coupling strength between the atoms is given by the 
coupling constant g. Backscattered photons are detected in 
four polarization channels: (a) right circular in, left circular 
out (with respect to a fixed observation direction; h \\ h); (b) 
linear in, orthogonal linear out {lin ± lin); (c) right circular 
in, right circular out [h L h); (d) linear in, parallel linear out 
{lin II lin). Note that in cases (c) and (d) CBS appears on 
the background of single scattering from independent atoms. 



k^. In the case of linear polarization, the quantization 
axis is conveniently chosen along the incident polariza- 
tion vector eo (perpendicular to ki,). For both choices, 
e has, in general, three nonzero projections on the unit 
vectors of the spherical basis, at finite angles between k 
and — k^. 

From now on, we will use an approximation which is 
legitimate at small scattering angles (k is very close to 
— ki), a common situation in CBS experiments. One 
is interested in signals at very small angles 9 ~ 1/ki 
around the backscattering direction (the mean free path 
i is the average distance between consecutive scatterers, 
and fc£ ^ 1 in dilute atomic gases) . Since the geometric 
change of e varies only with the cosine of the scattering 
angle ^ 1, we can take e constant, equal to its value 
at exact backscattering. 

With these conventions, Eq. H14|l can be specialized 
for the four polarization channels traditionally selected 
in CBS experiments (see Fig.|21l. 



1. Helicity preserving channel {h \\ ti) 

In the helicity preserving channel, the incident radia- 
tion is circularly polarized, and drives either the |1) — > |2) 
or the |1) |4) transition. The backscattered light 
is then observed in the orthogonal polarization channel 
(with the same helicity, since the propagation direction 
is reversed) and must be radiated by the |1) |4) or 
the |1) |2) transition, respectively. Both combinations 
are completely equivalent, and for the case El = e+i, 
£ = e_i, the total backscattered intensity reads 

I = E('^22>ss + Kl'^f2)s.e*-^°^ (15) 

a a#/3 
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2. lin _L lin channel 

In the lin _L lin channel, £l = gq. We further assume, 
without loss of generality, that the laser field is propa- 
gating along the x-axis, such that the detected photons 
are polarized along the y-axis. In the spherical basis, 
e = By = i(e_i + e+i)/\/2, leading to the expression 

I = ^((a'22)ss + (o-44)ss + (o-24)ss + (o-42)ss) 

a 

+ 5]e^'^-"-°^(Ki<),, + (a2"iaf2).. 

+ {<A2U + {aqAs). (16) 

3. Flipped helicity channel [h _L h) 

With incident polarization el ~ e+i as before, the 
h J- h channel corresponds to e = e+i, such that 

I = E«4)ss + E Kyidsse'""""^ . (17) 

4- lin \\ lin channel 

Finally, in the lin || lin channel, the incoming and 
outgoing photons are linearly polarized along the same 
axis, £l = £ = §0, and we obtain 

I = E('^33)ss + E KAse'"^-"-^ . (18) 

The above expressions systematically decompose into 
I 



two parts: Single atom, steady state dipole expecta- 
tion values express the intensities radiated by individ- 
ual atoms, while correlation functions of distinct atomic 
dipoles, multiplied by phases which depend on the rel- 
ative position of the atoms, account for the interference 
contribution. 

We now show how to evaluate these different steady 
state expectation values explicitly, before performing the 
ensemble average over the atomic positions, in Sec. IIIII 



C. Equations of motion for the atomic correlation 
functions 

The dynamics of the atomic dipole operators' expec- 
tation values as well as of the dipole-dipole correlators 
which enter (|15ll8f) is governed by the master equation 

m 

N N 

{Q) = J2{£aQ)+ E i^c^pQ), (19) 

where the Liouvillians Ca and Cafi generate the time 
evolution of an arbitrary atomic operator Q, for inde- 
pendent and interacting atoms, respectively, and (...) = 
Tr(. ..p(0)), with p(0) being the initial density opera- 
tor of the A'^-atoms-field system, indicates the quantum 
mechanical expectation value. 

In the co-rotating frame with respect to the driving 
field at frequency ujl, and after the standard electric 
dipole, rotating wave, and Born-Markov approximations, 
Ca and La/s read [l5| : 



-td[Bi ■ D„, Q] - -[n^iBl ■ el) + KiBa ■ el), Q] + 7 • [Q, D„] 



CapQ = Dl- T(.g,n).[Q,D;3] + [ 



Q] - T*(5,fi).D„ 
I 



(20) 
(21) 



where the radiative dipole-dipole interaction due to ex- 
change of photons between the atoms is described by the 
tensor T (5,6) = 7gA. This interaction has a certain 
strength depending on the distance between the atoms, 
via 

3f 



(22) 



with ko — luq/c, and on the life time of the excited atomic 
levels, through 7. 

The coupling constant \g\ <C 1 is small in the far-field 
{korap S> 1), where near-field interaction terms of or- 



der (kgrap) ^ and (fco^Q^) can be neglected (which, at 
higher atomic densities, could also be retained in our for- 
malism). The projector A = 1 — nn on the transverse 
plane defined by the unit vector fi along the connect- 
ing line between the atoms a and P is explicitly given 
through 



1 = -e_ie+i + eoGo - e+ie_i. 



(23a) 



e'-^sini? „ e-^^sini}^ , 

n = -= — e_i+cost/eo -j= — e+i (23b) 

V 2 V 2 

The angles {'&,((>), which fix the direction of the con- 
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necting vector between two atoms (with respect to the 
backscatterin g di rection), will have to be averaged over 
further down 

It should be kept in mind here that the master equa- 
tion treatment implies a trace over the modes of the 
free field, and that the Markov approximation implies 
some coarse graining on the time axis. Thus, 7 and g 
are the only remnants of the coupling to the electromag- 
netic vacuum, giving rise to some effective dynamics of 
the atomic operators, on time scales which are long with 
respect to the time scales of single absorption and emis- 
sion events from and into the electromagnetic reservoir. 
Only by expansion of the solutions of H20I21() in powers of 
g will we be able to distinguish multiple scattering con- 
tributions of increasing order, since, formally, all elastic 
and inelastic processes are lumped together in H2UI21|I by 
f2, 7, and g. This renders the master equation treat- 
ment somewhat less transparent or at least less intuitive 
as compared to the scattering theoretical approach [lTI |. 
but bears the qualitative improvement of yielding results 
which are valid for arbitrary saturation parameter s. 

Note that the Markovian master equation H19|l ignores 
retardation effects due to a finite photon propagation 
time between scatters. This approximation is justified 
when max(ra/3) cx ^ <C c/7 15]. In typical experiments 
with sharply defined, resonant optical dipole transitions 
P,0], both the condition of diluteness, k£ ^ 1, and of 
'instantaneous' propagation are very well satisfied. 

The equation (|19|l leads to a system of linear, coupled 
differential equations with constant coefficients for the 
atomic correlation functions. The algebra of the iV-atoms 
operators is spanned by tensor products of individual op- 
erators cr^; = \k)^ each acting on the iV-fold tensor 
product of the four dimensional Hilbert space in which 
the internal states of a single atom are represented. For 
the free evolution of a single atom {N — 1), Qa can be 
chosen in a complete orthonornial set of 16 basis opera- 
tors (see also Sect. Ill Dl below, for details). 

For our treatment of CBS, we need to include at least 
two-atoms operators in (|19|l . For N = 2, the number of 
equations in H19|) is 255 = 16^ — 1 (there is one constant 
of motion) . In matrix notation, the resulting equation of 
motion reads 

(Q) = (A -f V)(Q)-|-j, (24) 

where the elements of the vector (Q) are given by the 
expectation values of the complete orthonormal set of 
two-atom operators. The elements Anm, VJi™ and j„ of 
the matrices A, V, and of the vector j, are derived from 
the equation for the element {Qn) in 124|) : 

255 

{{Co, + Cp)Qn) = ^ A„„(Q™) -)-j„, (25) 

m— 1 

255 

((^Q/3 + ^f3a)Qn) — ^ VnmiQm)- (26) 
m— 1 

From the decomposition of H24|) into H25|l and H26|l it is 

apparent that the matrix A generates the evolution of 



uncoupled atoms, whereas V describes their interaction 
via the exchange of photons. 

D. Green's matrix and matrix A 

In order to solve (|24|l . we first need its represention in 
a suitable operator basis. Thereafter, we will derive a so- 
lution for independent atoms (that is, we will ignore the 
matrix V, which mediates the interatomic correlations) 
by a Laplace transform. The thus established relation 
between the Green's matrix for the non-interacting two- 
atom system and the matrix A will then serve as a basis 
for a systematic, perturbative treatment of the interact- 
ing case, at increasing order in the coupling constant g. 

Let us first consider the dynamics of a single four-level 
system. An expectation value of a single-atom operator 
Qa obeys the equation of motion 

(Qa) = {Co.Q»), (27) 

where the superoperator Ca is given by eq. H2U|) . The 
master equation H27|) describes resonance fluorescence of 
the laser-driven = — > Jg = 1 atomic dipole transi- 
tion. Qa belongs to the complete orthonormal set S"" of 
16 operators for the four-level system, 

12 operators 

(28) 

where 

1 = CTii + 0-22 + 0-33 + 0-44 , (29a) 
fJ-l cr22 - 0-33 + (J44 - (111 , (29b) 

A*2 = cr22 - 0-33 - CT44 + an , (29c) 
A*3 = 0-22 + 0-33 - (T44 - 0-11 • (29d) 

It is easy to check that for the elements of 5" the or- 
thonormality condition Tr [QnQm] = ^nm holds. In this 
representation, equation l|27|l turns into a linear matrix 
equation for the vector {Qa{t)), whose 16 elements are 
the quantum mechanical expectation values of the ele- 
ments of 5". Since the atomic levels' dynamics are un- 
coupled, except for the laser-driven transition, it can be 
solved analytically. The dynamics of the driven transi- 
tion is equivalent to the one of a two-level system. 

For two atoms, Q"^ S 3°" (g) S^, and the vector of the 
two-atoms correlation functions, 

(Q"/3(t)) = [(1" « l/5)/4, . . . , (aS ® af,)]^ , (30) 

consists of 256 elements. The evolution operator of un- 
coupled atoms reads 

e^°* = e^°* e^"*. (31) 

A Laplace transform dte~^*e^°* of H31|l gives the 
Green's function (or resolvent) Gafi{z) = {z — £o)~^ of 
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the Liouvillian Cq = + In the two-atom basis 
S°' (g) S^, Gapiz) has a matrix representation Gol(}{z). 
This matrix has the following block structure: 



Vz-ij : Go(z)/ 



(32) 



where vectors (zero vector) and j have 255 elements, 
and G[){z) is the truncated (255 x 255) Green's matrix. 
As seen from (j22J, the first column of matrix Gap{z) 
has a pole at z = 0. This pole appears because the first 
element of the vector (Q"'5(t)), (1" ® l'') = Trpo, is 
a constant of motion. All other elements of the vector 
(Q"''(t)) are time-dependent. The steady-state solution 
of the truncated vector, (Q(i)), which is obtained from 
(Q"''(t)) after exclusion of its first element, Trpo, is de- 
fined as (Q)ss — liniz—o -z(Q(^)), where (Q(z)) is the 
Laplacian image of the vector (Q(i)). This limit is eval- 
uated to give 



(Q)ss — Goj, 



(33) 



where Go = Gq (0) . Comparison of with the steady- 
state solution of Eq. H25|l. (Q)ss = — A^^j, now shows 
that Gn = -A^i. 



E. Perturbative restriction to low scattering orders 

The theoretical description of coherent backscattering 
is relatively simple in two opposite regimes: either in 
the diffusive regime of fully developed multiple scatter- 
ing in optically thick media, where long paths or high 
scattering orders yield the celebrated conical line-shape 
of the CBS signal, or, on the contrary, in the regime 
of scattering by optically thin media (or in the presence 
of suppression of interference), where only double scat- 
tering needs to be considered 0, 0|. Indeed, it is in 
the double-scattering regime where the first experimen- 
tal observation of a CBS reduction due to the non-linear 
saturation of atomic dipole transitions was reported 6]. 
We limit our present analytical and numerical analysis 
to this specific case. 

In our master equation framework, all information on 
multiple scattering processes is contained in the corre- 
lation functions of dipole-dipole interacting atoms. The 
double scattering contribution to the scattered light in- 
tensity from a given pair of atoms, resulting from the 
exchange of two photons between the atoms, is obtained 
by perturbative expansion of the respective correlation 
functions to second order in the dipole-dipole coupling 
constant \g\. This contribution depends only on the ob- 
servables related to the two selected atoms, and not on 
those of the rest of the atoms in the cloud, since correc- 
tions to the mean intensity due to the latter would be 
of higher order in \g\ {\g\^ , etc.). Therefore, to find the 
double scattering contribution, we will solve the master 



equation H19() for two fixed atoms a = 1,2. Subsequently, 
we have to add up all double scattering contributions re- 
sulting from the atoms located at random positions. In 
other words, we need to perform appropriate disorder 
averages of the solution for two fixed atoms. 

This setup defines the simplest possible model describ- 
ing CBS. Indeed, double scattering is the lowest order 
process which gives rise to distinct scattering amplitudes 
that can interfere constructively. Despite its simplicity, 
this model allows for a qualitative assessment of the im- 
pact of nonlinear scattering processes on the CBS signal, 
whereas propagation effects in the bulk of the scattering 
medium are beyond its reach. It also needs to be consid- 
ered that nonlinear scattering processes are induced by 
high laser intensities, at which atoms are rapidly acceler- 
ated out of resonance. Within our model, we neglect this 
acceleration, by focusing exclusively on the coupling of 
photons to the internal atomic degrees of freedom. Such 
an approximation is justified, since the mechanical action 
of light on atoms can be experimentally compensated by 
shortening the CBS probe duration, as realized, e.g., in 
i- 

For our perturbative solution of Eq. H24|) , we take ad- 
vantage of the small parameter g in \'2,1\ (defined in 
Eq. 1)22(1 '). and expand in a power series of V. The nth 
order 



(Q)N = (GoV)"Goj 



(34) 



of the stationary solution of Eq. (|24|l gives the two-atoms 
correlation functions resulting from n exchanged pho- 
tons, and generally includes also recurrent scattering (a 
photon visits the same atom several times). We recall 
that, in the regime of elastic scattering from dilute sam- 
ples of resonant scatterers, higher scattering orders can 
be accounted for by considering more scatterers, while 
recurrent scattering is irrelevant [l9| . In contrast, as 
one proceeds to the strong scattering regime, with denser 
clouds of resonant scatterers, the effect of recurrent scat- 
tering manifests by a gradual reduction of the enhance- 
ment factor as compared to its maximum value 2.0 |2Clj| . 
We wiU see in Sect. IIV C II below that Eq. gener- 
ally including recurrent scattering contributions is fully 
compatible with neglecting these in the linear regime. 

Total intensities for the double scattering contribution 
are given by the third term (proportional to \g\'^) of the 
above expansion of the correlation functions which enter 
Eqs. (|15I18|) . with the general structure 



(Q)[21 = GoVGoVGoj. 



(35) 



Note that the matrix V depends on A , though neither 
on El nor on e. Indeed, information about the laser po- 
larization is carried by the Liouvillian l|2U|) , which governs 
the evolution of independent atoms. As for e, it does not 
appear in the equations of motion for the atomic cor- 
relation functions at all - but it defines which elements 

[2] 

of the vector (Q)ss contribute to the observed signal, as 
evident from Eq. (|14|l . 
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F. Elastic component of double scattering 

The total backscattered intensity H14() has a spectral 
distribution that contains an elastic and, beyond the 
weak field limit, also an inelastic component. The de- 
tected intensity is the correlation function of the source 
field amplitudes radiated by the atomic dipoles. 

Its elastic component stems from the classically radi- 
ating dipoles, i.e., from the nonfluctuating factorized av- 
erages (f^j) 0. Hence, the elastic intensity is given 
by the product of the expectation values of the atomic 
dipoles, 

N 

^ (£.Di),,(£*.D0),,e*-°''. (36) 

a,/3=l 

In the helicity preserving channel h \\ h, this reads 

- E K'^2i>ss|^ + E«)-(^f2)s.e^''-^°^. (37) 

Analogous expessions for the elastic component can be 
derived from Eqs. H16I18|I in the other channels. 

The double scattering contribution /|' to the elastic in- 
tensity, which we are interested in, and which is propor- 
tional to |f;p, is obtained from the corresponding power 
series expansion, Eqs. (|34I35|I . of the individual factors 
entering the above expression, in the coupling constant. 
In this expansion, there emerge symmetric and asymmet- 
ric combinations, like (o-2i)is' (CTi2)is , and (o-2i)is' (a'?2)£' , 
respectively. To ease the physical interpretation of these 
various terms, remember that, by virtue of (|34|l . the su- 
perscripts '[0]', '[!]', and '[2]' signal the scattering of a 
photon from one single atom a ('[0]'), subsequently from 
atom a and then from atom (3 ('[!]'); and the rescattering 
of the same photon from atom a, after a first encounter 
with a and subsequent scattering from (3 ('[2]'). 

In the h \\ h and lin _L lin channel, only symmetric 
combinations contribute to the elastic component of the 
CBS intensity, since the lowest order expectation value 
of a single atom's coherence in the analyzed transition 
(like ( '^i2)ss in the helicity preserving channel) must van- 
ish: at lowest order, the atom is not coupled to the other 
atom, nor is the transition directly driven by the injected 
laser. By the same argument, also the single scattering 
intensities from non-interacting atoms (arising from the 
first sum on the right hand side of Eq. H37|l ) are projected 
out, and an interference signal from purely multiple scat- 
tering sequences is measured in these channels (unless the 
atoms have a degenerate ground state and can undergo 
transitions between different Zeeman-sublevels |2l|'). 

In the h J- h and lin \\ lin channels, both, symmetric 
and asymmetric combinations of products of the dipole 
averages contribute. 

Altogether, at second order in the coupling constant, 
the elastic backscattering intensity from two fixed atoms 
is obtained from the evaluation of the zeroth-, first-, 
and second-order stationary solutions (Q)is', (Q)is', and 



(Q)ss of Eq. if^ . The backscattering intensity from 
a cloud of randomly located atoms is finally derived 
through an appropriate disorder average. 



III. DISORDER AVERAGING 

Coherent backscattering is such a surprising effect be- 
cause it survives the ensemble average over random posi- 
tions of the scatterers, which destroys all other speckle- 
like interferences. The precise procedure of disorder av- 
eraging is important as soon as one is interested in the 
exact shape and angular width of the CBS cone. Here, 
we rather focus on the impact of the atomic saturation 
on the maximum CBS intensity, in the exact backscat- 
tering direction. Since saturation effects are independent 
of the precise averaging prescription, we choose a proce- 
dure as simple as possible: an (i) isotropic average of the 
relative orientation n of the atoms over the unit sphere 
is followed by (ii) an average of the inter-atomic distance 
ri2 over an interval of the order of the laser wavelength, 
around their typical distance i: 

^ l-e+27r/kL r 
(...)conf = / dri2 dnn.... (38) 

47r Ji>-27T/ki^ J 

After the evaluation of H38|l . the final expression for 
the backscattering intensity has a general structure which 
decomposes into 'ladder' and 'crossed' contributions, re- 
spectively. The ladder terms collect the intensities scat- 
tered by individual atomic dipoles (in a diagrammatic 
representation, they arise from the summation of co- 
propagating amplitudes 0] along a sequence of scatter- 
ers), in an incoherent sum, whereas the crossed terms 
stem from the interference of amplitudes radiated by dis- 
tinct dipoles (counter-propagating amplitudes in a dia- 
grammatic picture) , and are garnished by the associated 
phases. Through this phase factor, the interference part 
depends on the angle 9 between the wavevector k of the 
final photon and the backscattering direction — k^. As 
an example, the ladder and crossed terms in the h \\ h 
channel read, by virtue of Eq. H15|l . 

LT - (('^22)^1 + (^22>''l)..>conf, (39) 

CnO) = 2Re((a2\a?2)Se^'^-"--)eonf, (40) 
respectively. 

The enhancement factor a, which is the figure of merit 
for the quantification of CBS, is given by 

^2 

An analogous expression for the elastic CBS component 
follows from Eq. H37|) . with ladder and crossed terms Lf 
and Cfie). 
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FIG. 3: (Color online) Saturation dependence of the total 
double scattering contribution /2°' to the CBS signal in the 
/i II /i channel, decomposed in its ladder and interference parts 
1/2°' and C2°', according to Eqs. 146I47II . at exact resonance 
(5 = 0. At finite saturation s, the interference term C2°' drops 
below the ladder contribution Z/2°' i indicating a loss of coher- 
ence. At large s > 1, the total double scattering intensity 
must decrease with s, since the scattering cross section of the 
emitting atom drops as . 



i?i(s), i?2(s), and P{s) are polynomial expressions in the 
saturation parameter s, 



(6912s + 3168s2 



9 

+2645^ - 
i (1152s 



n , (44a) 
528s2 + U2s^ + 7s*) , (44b) 
(l + s)2(12 + s)(32 + 20s + s^), (44c) 



and 



"9' 7 



{q,q'^ ±1,0). 



(45) 



The configuration average over H43|l and 1)42(1 . defined 
in Eq. (|38|l . leads to the final result 



r tot 
^2 



\~9\'Rii-s 



(4 + s)P(s) 
15P(s) 



) /2 (Hg)^ 
s)U5 35 /' 



(46) 
(47) 



IV. RESULTS 



We now proceed to evaluate the general expressions 
derived above, for the four typical polarization channels 
analyzed in the laboratory: h \\ h, lin _L lin, h ± h, and 
lin II lin. Analytical results for vanishing detuning 5 will 
be complemented by some numerical results for 6 ^ Q. 



with g = g\r^ij=i (see Eq. (|22|) '). The scattering angle 
6 = 2arcsin{|k + kL|/2fcL} ^ 1 with respect to the 
backscattering direction was assumed to be sufficiently 
small herein. A power series expansion of Eqs. ((46I47|I 
to second order in s reproduces the diagramatically ob- 
tained result of llj, 



^tot 



cx s ■ 



L*2°' oc s 



(48) 



for the total double scattering intensity, at 6 = 0. 

The behavior of C^°*(0) shows that the dou- 

ble scattering intensity = L^°* -I- C^°*(0) behaves 
markedly different from that of an isolated atom. While 
the radiated intensity from an isolated atom. 



h \\ h channel 



j[0] 



1 



(49) 



1. Total intensity, at zero detuning 



In Eqs. (ES)) and (|2Z|) we assumed that the |1) |4) 
transition is laser-driven. Photons with preserved helicity 
originate from the |2) |1) transitions, and £l = e+i, 
£ = e_i. The total double scattering intensity for two 
fixed atoms reads, by virtue of Eqs. H15() and H35() . at 
second order in g and vanishing detuning 6 = 0: 



2Re {{al.al,: 



[2]„Jk-ri2 



Rii-s) 



2 \[2] 



(4-f s)P(s) 
xcos{(k + ki)-ri2}, (42) 

2| 
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IA+t,.iP^. (43) 



Pis) 



grows monotonically with s until it finally saturates [3], 
the double scattering intensity exhibits a maximum at 
s ~ 0.7 (see Fig. 13), followed by gradual decrease oc s^^ 
for large s. Also this is a simple consequence of the 
saturation behaviour of an isolated atom described by 
(|49|l : At high injected laser intensities, the atom that 
emits the final photon has a total scattering cross sec- 
tion that asymptotically decays like I^^^/Il oc s~-^ (see 
Eqs. Hllll|l 'l. and is consequently less likely to scatter 
photons coming from the other atom. 

The enhancement factor a{s), Eq. (|41|l . deduced from 
Eqs. (|46I47|I reads 



a{s) = 1 



Riis) 



(4 + s)i?2(s) 



(50) 



and a(0) = 2.0 in the weak field limit, as expected. The 
dependence of a on the saturation parameter is shown 
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FIG. 4: Numerical enhancement factor a in the hehcity pre- 
serving channel h \\ h, versus saturation s of the atomic 
transition, for different detunings 5 from resonant driving. 
Larger detuning leads to a faster loss of the CBS contrast. 
At small s, one recovers the perturbative prediction fil| lin- 
ear in s (straight dashed line). At large s, a saturates at 
a 5-dependent value as^ > 1, due to the constructive self- 
interference of inelastically scattered photons. Remarkably, a 
also passes through a minimum at s ~ 0.5, for very large de- 
tuning S = 2O7. This indicates destructive interference, and 
the physical cause of this observation remains to be identified. 



in Fig. 2] As above for the individual cross and lad- 
der terms, we again obtain perfect agreement with the 
linear decay predicted by the scattering theoretical re- 
sult a ~ 2 — s/4 11], in the limit of small s. When s 
increases further, a monotonically drops to an asymp- 
totic value lims^oo Q^(s) — Q^oo = 23/21 which is strictly 
larger than unity, implying a nonvanishing residual CBS 
contrast in the limit of large injected intensities. As we 
shall see further down in Sec. IIV A 31 this residual con- 
structive interference effect is exclusively due to the (self- 
)interference of inelastically scattered photons. 



2. Finite detuning 

It is in general no more possible to obtain explicit ex- 
pressions for the Green's matrix Go (tantamount of in- 
verting A), in the case of nonvanishing detuning 5^0. 
However, this can always be done numerically, and Fig.0| 
compares the enhancement factor at resonance to the one 
for three different nonvanishing values of 6. In qualita- 
tive agreement with the experiment 0, a decays faster 
for larger detuning, as s is increased from zero. But not 
only does the enhancement factor exhibit a steeper (ini- 
tial) decrease with S: it also reveals destructive interfer- 
ence (a < 1) for large detuning 6 — 2O7, at s ~ 0.5. This 
corresponds to a large Rabi frequency fl ~ 2O7. 

To gain some insight on whether such destructive in- 
terference is generic for large Rabi frequencies and large 
detunings, we monitor enhancement factor vs. detun- 
ing, for two fixed, large values of the Rabi frequency, 
as displayed in Fig. [S] For a given value of il, the en- 
hancement factor decreases as a function of \S\, from its 




20 30 40 50 60 
151 (in units of y) 



90 



FIG. 5: (Color online) Numerical enhancement factor a (nu- 
merical solution) in the helicity preserving channel, versus the 
laser detuning S, for two values of the driving Rabi frequency 
Q. The decrease of a below unity for = 2O7, in the vicinity 
of \S\ = 2O7, correlates with the minimum displayed by a(s) 
in the corresponding plot in Fig.0J at s ~ 0.5. 



maximum at (5 = to its f2-dependent minimum value at 
\S\ ~ Q. For very large detunings, as saturates at a level 
1 < lim|i-|^o2 = as^ < aoo, indicating (i) constructive 
(self-)interference of far-detuned photons, and (ii) similar 
behavior of the ladder and crossed terms, asymptotically 
in \S\. The asymptotic value as^ is the lower the larger 
fl. Furthermore, as a direct counterpart of the destruc- 
tive interference observed in Fig. ^ a drops below unity 
for fl — 2O7, in a finite range of \S\. 

Note that a similar effect was predicted in for 
linear double scattering from atoms with Zeeman-shifted 
hyperfine ground levels. In our case, the onset of de- 
structive interference at (5 ~ I57 and sq ~ il^/27^ = 
200 occurs approximately at saturation s ~ 0.9. A 
physical interpretation of this interference-induced anti- 
enhancement of CBS, which we tentatively attribute to 
an AC-Stark shift of the laser-driven atomic sublevels, 
will require a closer inspection of the total stationary in- 
tensity, and is refered to a separate contribution. 



3. Elastic component at finite detuning 

To see that the residual contrast observed in Fig. ^ for 
large saturation parameters stems from inelastically scat- 
tered photons, we now derive expressions for the elastic 
ladder and crossed contributions to the double scattering 
CBS signal. To do so, we extract the elastic contribution 
/j' to the total double scattering intensity from an ex- 
pansion of Eq. (|57|l to second order in \g\, as prescribed 
by ^MB- 



ik-ri2 



(51) 
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With the evolution equations 



{-7 + iS){a1^)--n^{a^^) 



2J 



-27K2> - ifll 



'12/ 



(52a) 



(52b) 



which can be derived from Eqs. H19I21|I . an analytic ex- 
pression for the steady state mean value (cr"2)ss' and 
thus for /|' is obtained by the following argument (note 
that this remains valid also for finite detuning S): As 
long as we content ourselves with a lowest order treat- 
ment of multiple scattering effects, only factorized zeroth- 
order correlation functions (Qf )£' (Q^ )s? for indepen- 
dent atoms a 7^ /3 contribute to the sums on the rhs 
of Eqs. jH2al52bll . since the coefficients are of order 
\g\, expressing the dipole-dipole interaction between dis- 
tinct atoms. Such products vanish except if their factors 
involve only the driven levels |1)^ or |4)^. Consequently, 
the summations in Eqs. (|52a|l and l|52bll . which extend 
over different subsets of the two-atom correlation func- 
tions, can be condensed according to 



19 



x(af4)M(a?i)M, 



(53) 



i,3 



75* ■^+1,+! 



x(af4)M«i)i°l, (54) 

respectively. Substitution thereof into (|52aj52b|l (with 
the Ihs of (|52al52b|l equal to zero, and {(t"^) — {a'^^)*), 
together with the known solutions of the optical Bloch 
equations for a single two- level atom 



« \[0] ^ iSl 

'14/ss 



i5) 



leads to 



'l2/ss 



1 + s' 



Ki) 



[0] 



2(1 + s) 



«75 



2(7-i5)2(l 
If we now rewrite Eq. H51|) as 



-Trc 



with Tdi, 



cjlof^e-'"^-^^ andT,, 



(55) 



(56) 



(57) 



(a?,)We-'''-^ the 



elastic component of the double scattering intensity ap- 
pears, with (|ll|l . as the square modulus of a sum of the 
'direct' and 'reversed' scattering amplitudes 



Tdir ~ 



«75 



2(7-iJ)2(l + s)2 
2(7-i,5)2(l-hs)2' 



(58) 
(59) 



Direct and reverse amplitude are symmetric under the 
interchange k ^ — k^, and thus satisfy the condition of 
reciprocity Consequently, the elastically scattered 

photons remain strictly coherent, for any s, and show per- 
fect CBS contrast with a{s) = 2, as immediately spelled 
out by the explicit expression for the elastic ladder and 
crossed terms, at arbitrary detuning and Rabi frequency: 



rcl 
^2 



2 1.9 



;;|2 



15 1 + (,5/7)2 (1 + sY 



(60) 



While perfectly coherent even for large s, the elastic 
contribution to the total double scattering intensity de- 
creases as s~^, by virtue of Eq. after passing 
through a maximum at s = 1/3. In contrast, the to- 
tal signal fades away like s~^, according to Eqs. H46I47|I . 
Consequently, in the large s limit, the CBS signal is com- 
pletely dominated by the inelastic scattering component, 
and the residual CBS contrast a^a observed in Sec. lIV AH 
above is due to the selfinterference of inelastically scat- 
tered photons, which are incoherent with respect to the 
injected laser radiation. The visibility of this residual in- 
terference signal is limited by the amount of which-way 
information communicated to the environment, during 
the multiple scattering process |23| • 

Finally, let us note that equation (|60|l allows for a 
transparent interpretation, since it can be factorized into 



(i) the elastic intensity 

jol[0] ^ 



(1 + S)2 

scattered by the first strongly driven atom, 
(m) the total scattering cross section 



(61) 



(62) 



(l + (V7)^)(l + s) 
of the second atom, and 
(Hi) the relative weight /ci[o]/jtot[o] ^ cr'^Vo''"* equal to 



1 



of elastic processes 



^2 



(63) 



B. lin ± lin channel 

Up to a constant factor 1/2, the results for the lin _L 
lin channel turn out to be the same as for the h \\ h 
channel, at exact backscattering. For a given value of s, 
the ladder and crossed terms are two times smaller than 
in the /i || /i channel. Since, however, Eq. (|16|l for the 
intensity in the lin _L lin channel is manifestly different 
from Eq. H15|) . the h \\ h result, a short discussion of this 
observation is in order. 
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In both cases, the CBS intensity is observed in the 
polarization channel orthogonal to the excitation chan- 
nel. In the h \\ h channel, the orthogonal channel 
is defined by one dipolc transition |2) ^ |1). In the 
lin _L lin channel, the orthogonal channel is defined by 
two atomic transitions, |2) |1) and |4) |1). Yet, 
by introducing a superposition state |e) = |2) + |4), we 
can rewrite expression (|16|l in a way which is formally 
equivalent to (|15|l . However, the geometric weight of 
the resulting ladder and crossed intensities is given by 

( Ao -1+ Ao,+i)( A+i,o+ A_i^o)/2. Only two of these 
four terms, A o.~i A _i.o/2 and A o,+i A +i.o/2, survive 
the configuration average, what leads to results that are 
two times smaller than in the h \\ h channel. 



C. h A- h channel 

We shall now consider detected photons which have 
the same polarization as the incident ones. Hence, as al- 
ready briefly discussed in Sec. IIIFI single scattering as 
well as double scattering events will contribute to the 
detected signal, with the latter only a small correction 
to the former. Correspondigly, an experimental detec- 
tion of the double scattering contribution alone is ex- 
cluded. Nonetheless, it is instructive to consider this sce- 
nario in the regime of a nonlinear atomic response to the 
injected radiation because it allows to identify the role 
of recurrent scattering where one photon rescatters from 
the same atom, after visiting the other (this process is 
second order in the coupling constant \g\). 



1. Total intensity 

To obtain explicit expressions for the total scattered 
intensity, we proceed stepwise and first expand the single 
atom contribution to the intensity, in Eq. H17|) , to second 
order in g. We obtain 



c«=l 



where 



i^i(s)|A 



-1-1, 



(64) 



+i^2(s)(|A+i,oP + |A+i,+in 
-t- terms cx Re [g"^] , \g\^ cos{2kL ■ ri2} , 



36s + - 27 - 19s-^ - 
12(l + s)5(3 + s) 
4s2(288-M32.s + 23s2 + s3 
3(1 -hs)P(s) 



(65) 
(66) 



with P{s) from Eq. The terms in the last line of 
oscillate rapidly on the typical scale £ of the interatomic 
separation (recall Eq. if^ . and k^rap ^ 1; this is, also 
the dependence of g on ri2 is to be taken into account 



here), and average out under the integral over ri2 in 138|l . 
Thus, they will be dropped hereafter, whereas terms oc 
\g\^ = 9/4fco?'i2 vary smoothly with ri2, and will be kept. 

An analogous expansion of the interference terms in 
(|17|l yields the expression 



1.91^^^3(5)1 A+i,_i|^ (67) 
X cos{(k -I- ki) • Y12} 
+ terms cx Re [g"^] , 
|gpcos{(k-kL) •ri2}, 
|gpcos{(3ki~k) •ri2}. 



where 



324s + 540s^ + 450s'* + 2l%s'^ + BSs^ 



29s^ 



36(l-hs)6(3-hs)2 



(68) 

The three last lines of Eq. H67(l are irrelevant for our 
subsequent treatment, for exactly the same reason as the 
corresponding terms in Eq. H64() . 

We now have a closer look at the geometric factors in 
these equations, which allow the identification of the un- 
derlying elementary scattering processes. The geometric 
weight of the contribution proportional to F2(s) in (|64|l is 
given by | A +i_oP + | A ^i^+ip. It describes the couphng 
of a photon from the |4)^ |1)^ transition to either the 
|3)/3 |l)/3, or to the \2) ^ — > |1)^ transition, respec- 
tively, then back to the |4)^ |1)^ transition, and only 
then to a detector. This is recurrent scattering. (Note 
that non-recurrent transitions, e.g., from |4)^ |1)^ to 
|3)^ — > or to \2)p cannot give rise to a de- 

tected photon in the h J- h channel.) 

More precisely, the single scattering contribution with 
the weight F2{s) originates from the interference between 
single scattering from independent atoms and (recur- 
rent) triple scattering, in which a photon is subsequently 
scattered by atom a, then by atom (3, and by atom a 
again (recall our discussion in Sec. IIIFI and also see 
Eqs. (|72I73I) in our subsequent discussion of the elastic 
contribution to the /i _L /i channel). Since the rate of 
single and recurrent scattering is equally limited by the 
number of photons incident on the atom, cx s, it follows 
that F2{s) cx s^ for s — > 0. This is consistent with a ba- 
sic postulate of multiple scattering theory jl^ , according 
to which recurrent scattering is irrelevant in the linear 
regime of weak saturation. It is also clear why there is 
no recurrent scattering in || /i and lin _L lin channels: 
Indeed, single scattering is essential for recurrent scat- 
tering to show up at order but is filtered out in the 
orthogonal polarization channels. 

Now consider those terms in Eqs. (|64|l and H67I) with 
angular part |A_(.i^_ip. Along with the processes in 
which atoms exchange photons, these equally much rep- 
resent the two recurrent scattering sequences 
|4)/3 ~* \^)a- The weight of these contributions is not the 
same as for |4)„ ^ |3)^ ^ |4)„ and |4)„ |2)^ |4)„ 
transitions, since the presence of the driving field in the 
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FIG. 6: (Color online) Total crossed (Eqs. (jTOJ) and ladder 
(Eq. 1691 ') terms, together with the double scattering intensity 
Jt°' = L^°* + C2°'(0), in the h ± h channel, as functions of 
the saturation s, at exact resonance (5 = 0. 



|1) |4) transition definitely destroys the symmetry of 
the excited state sublevels. 

As regards the sign of the various scattering contri- 
butions in Eqs. H64I67() . note that the weight of the in- 
terference part given by Eq. (|68|l is positive for all s, 
while Fi{s) and ^2(5)1 Eqs. I|()5lt)()l) . are nonpositive, the 
function ^2(5) being strictly negative for s ^ Q. Hence, 
recurrent scattering is a small negative correction to the 
overall positive single-atom contribution to the total scat- 
tering signal. 

The total crossed and ladder contributions to the dou- 
ble scattering intensity are once again obtained after a 
final configuration average of Eqs. (|64|l and (|67|) : 



T tot 

^2 



l-\g\^(7F,{s) + 3F,{s))., 



(69) 
(70) 



The result is plotted in Fig. El where L^'^ turns nega- 
tive in a finite interval of s. In the limit of large s, the 
ladder term approaches zero from below, as 1/s, whereas 
the crossed terms decreases towards zero, with the same 
rate 1/s. Note that the negativity of the ladder term 
is compensated for by the single scattering contribution 
Li which cannot be separated from the double scatter- 
ing contribution, in the h J- h channel. Furthermore, 
the single scattering contribution does not decrease with 
growing s but rather saturates, so that for very large 
saturation parameters we can simply ignore the double 
scattering contribution. For all s, the CBS enhancement 
factor reads, after inclusion of Li: 



a = 1 



_C*°*(0) 



(71) 



In a real medium, the relative weight of single and 
double scattering depends on the optical thickness. How- 
ever, our simple model cannot correctly account for this 



effect; hence, we cannot assess here whether the negativ- 
ity of the double scattering ladder term has observable 
consequences in laboratory experiments. 



2. Elastic component 

Let us finally extract the elastic component of the CBS 
intensiy in the h J- h channel. According to (|17|l and 
with l|38|) . the elastic ladder and crossed terms are 
given by 



Cf{9) 



+ (-ll>i"sVi4>?l 

2 Re (e^''-- {{al)^^J{aW 



(72) 



-li)ff('^i4>i°l) , 

/ conr 



(73) 



\[0]/^2 W2] 



/ I \[2]/^2 x[0] 
\^4l/ss \''^14/ss 



Due to the factorization (|36|l of the classically radiat- 
ing dipoles, symmetric and asymmetric scattering con- 
tributions are directly born out: Products of first order 
(labeled by '[1]') contributions in g represent double scat- 
tering of photons subsequently at atoms 1 and 2, in direct 
and reversed order, whereas products of zero and second 
order ('[0]' and '[2]', respectively) express indistinguish- 
able single scattering and recurrent scattering amplitudes 
upon either one of the atoms. Explicitly, the above ex- 
pressions have the following geometric weights. 



rcl 

^2 



|2\ 



conf 



3r(|A+i,o| 
2F2(s) 



S + S''^ 

(1 + sY 



+ IA+1.+1P) 



Fi(s) 



conf 



l + s 

|2 



(74) 



Cf{6) = (|.9nA+i,_i|2cos{(k + ki).ri2} 



+ F5is) 



conf 



(1 



(75) 



where all terms which vanish under the configuration av- 
erage have already been dropped (see also the discussion 
of Eqs. (|64I67(I above), and 



F,{s) 



-36s2 - 39s3 - Us* + s' 

3(1 s)6(3 -f s) 
-72s2 - 51s3 - s'^ + 3s5 



(76) 
(77) 



Upon evaluation of the configuration average, we ob- 
tain the final result 



Li' - ir^ 

7 



15 



+ F5{s) 



6F2(s) 



7 1 



(1-t-s)^ 



,(78) 
(79) 
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Expressions (|78|l and H79|l imply that, in the small-s hmit, 
the elastic ladder and crossed terms coincide, as for or- 
thogonal polarization channels (see Eq. H6Q() ^. However, 
this equipartition does not prevail here beyond the lin- 
ear regime, due to reciprocity violating processes pjj 
which lead to a deviation of Lf from Cf{0), already at 
quadratic order in s. The violation of reciprocity was 
originally demonstrated in ^IL] for scalar atoms, and it 
will be demonstrated below in Sect. lIVEl that the scalar 
results immediately follow from our present results for 
parallel excitation/detection polarization channels, when 
ignoring those electronic sublevels that mediate recurrent 
scattering. However, beyond those reciprocity-violating 
processes already implicit in the scalar treatment, other 
processes specifically due to the vector character of the 
injected radiation field lead to additional deviations (ex- 
pressed by the term cx -^2(5) in t^78\i ). 

Furthermore, note that the elastic ladder and crossed 
intensities in the h J- h channel asymptotically behave 
like oc 1/s^ and oc 1/s, respectively. Also the total dou- 
ble scattering intensity in this channel is characterized by 
an asymptotic decrease oc 1/s (see Sec. IIVCTT)) . Conse- 
quently, unlike the h \\ h channel, where double scatter- 
ing becomes purely inelastic for large s, here both, elastic 
and inelastic photons, are present for large s. 



D. lin \\ lin channel 

The results for the lin || lin channel are the same as 
for the h J- h channel, modulo the following substitution 
of the geometric weights: 



'0,0| 



-1,01 



(80) 



This entails slightly different final expressions for the lad- 
der and crossed terms. For the total intensities, we obtain 



T tot 

^2 



1 



-\~g\^[8F,{s) + 2F,i, 

1 o 



CnO) = ye\~9\'F3{s) 
and the elastic result reads 



CfiO) 



l~|2 o 

1^1 15 



s + s^ F2{s) 
(1-1-5)6 + 2(1 + s) 



(81) 
(82) 

+ F4(s)] ,(83) 
(84) 



-(1-1-5)6 

Eqs. and differ from CTTHll and fT^m 

only through numerical coefficients. Therefore, all our 
above conclusions for the elastic and inelastic compo- 
nents of double scattering in the h J- h channel also apply 
for the present lin \\ lin case. 



E. Scattering of scalar photons on a two-level atom 

To conclude, let us briefly consider the model scenario 
of scalar photons scattering on a two-level atom - a wide- 



spread setting in typical quantum optical model calcula- 
tions, which neglects the important role of the polariza- 
tion degree of freedom in the presently discussed quan- 
tum transport problem. The scalar case is easily deduced 
from the above results for the /i _L ft, or lin \\ lin chan- 
nels, by simply setting equal to zero those contributions 
with the geometric weight | A+i^oP + I ~ since 

a two-level atom does not offer the required atomic tran- 
sitions. 

The total ladder and crossed terms then read, by virtue 

of Eqs. dMEZElZni), 



(85) 
(86) 



with the quadratic expansion 

L'2°* cx s - , C^°t(0)cxs-5s2, (87) 

in precise agreement with the result of [TH . 

Correspondingly, the elastic contributions, 
Eqs. (|78I79|I . reduce to 



Lf = 1^1 ^3 



a + s) 

s 

O + s) 



6 + F4is) 

,+F5{s) 



(89) 



which once again reproduces the second order expression 



C^'(O) oc s - 8s2 



(90) 



of Til. 

Thus, the scalar model correctly predicts the maximum 
enhancement factor a = 2, in the elastic scattering limit 
s ^ 0. Furthermore, it correctly describes those CBS 
contributions in the parallel excitation/detection chan- 
nels which originate from the laser driven transitions, for 
arbitrary s. However, the scalar model in general leads 
to incorrect results, since it ignores contributions from 
those sublevels of the degenerate excited state which are 
not driven by the laser, yet mediate recurrent scattering, 
as we have seen in Sec. IIV C II 



V. SUMMARY 

In summary, we have given detailed account of the 
master equation treatment of coherent backscattering of 
light from a disordered sample of cold atoms with a non- 
degenerate electronic ground state. This approach ex- 
tracts all physical observables from the steady state ex- 
pectation values of atomic dipole operators. Further- 
more, our treatment incorporates, to the best of our 
knowledge for the first time, the effect of interatomic 
dipole-dipole interactions for distant atoms. 
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In particular, the formalism allows to treat arbitrary 
pump intensities which possibly saturate the relevant 
atomic transitions, thus leading to inelastic scattering 
events (when more than one photon is incident on the 
scattering atom - on the spectral level, this entails the 
emergence of the famous MoUow triplet, in a single 
atom's fluorescence). The price to pay is a rapidly in- 
creasing dimension of the Hilbert space spanned by the 
many atoms' degrees of freedom. This limited our present 
treatment to two atomic scatterers, which is the mini- 
mum number of constituents to observe the CBS effect. 
Nonetheless, this approach allowed us to show that a 
small residual CBS signal survives even in the limit of 
purely inelastic scattering, due to the self-interference of 
inelastically scattered photons. 



Furthermore, we have seen that recurrent scattering 
leads to a reduction of the total double scattering signal, 
due to a destructive interference between single and triple 
scattering events upon the same atomic scatterer. 

Another advantage of the master equation treatment 
presented here, so far unexplored, is the immediate avail- 
ability of the CBS spectrum through a Fourier transform 
of suitable atomic dipole correlation functions, as well as 
of the associated photocount statistics. Whether CBS 
has an unambiguous signature in the spectrum and/or in 
the photocurrent remains hitherto an open question, but 
is getting in reach for state of the art experiments. 

It is a pleasure to acknowledge entertaining and en- 
lightening discussions with Dominique Delande, Benoit 
Gremaud, Christian Miniatura, and Thomas Wellens. 
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